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Abstract. By using a result from the numerical algebraic geometry package Bertini we 
show that (up to high numerical accuracy) a specific set of degree 6 and degree 9 polynomials 
cut out the secant variety 04 (P 2 x P 2 x P 3 ). This, combined with an argument provided 
by Landsberg and Manivel (whose proof was corrected by Friedland), implies set-theoretic 
defining equations in degrees 5, 6 and 9 for a much larger set of secant varieties, including 
0-4 (P 3 x P 3 x P 3 ) which is of particular interest in light of the salmon prize offered by E. 
\ Allman for the ideal-theoretic defining equations. 

<N . 

1 — 1! 1. Introduction 

O 

In 2007, E. Allman offered a prize of Alaskan salmon to anyone who finds the defining 
^ ■ ideal of the following secant variety 

: a 4 (p 3 x p 3 x p 3 ) , 

[A1110] . Recall that if A, B, C are vector spaces, then the Segre product is defined by the 
following embedding into the tensor product 

Seg: FA x FB x PC ->■ F(A ® B <g> C) 

([a],[6],[c]) ^ [o®6®c]. 

Further recall that if X C P^ is a variety, then the /c-secant variety of X denoted cr fe (X) C 
! P , is the Zariski closure of all points on secant P fe_1 's to X. For simplicity, we will drop 
the reference to the Segre embedding and write at{FA x FB x PC) for the secant variety 
to the Segre product. Secant varieties have been studied classically, but we have a renewed 
interest in their study because of the salmon prize and other related recent works on the 
subject (see |LM081lCl?U0c^lLan081[AR^[FriT0llLW071[S^09] ). 

Allman's ideal-theoretic question is still open. Our main result is Theorem 13.101 in which 
we give a geometric argument (relying on results of Landsberg and Manivel |LM08j Corol- 
lary 5.6] and the recent correction of the proof by Friedland in |Fril0j ) combined with a 
calculation using numerical algebraic geometry to show that up to high numerical accuracy, 
04 (P 3 x P 3 x P 3 ) is cut out set-theoretically by 1728 equations in degree 5, 1000 equations 
in degree 6 and 8000 equations in degree 9. Even though these dimensions are large, we show 
that in each degree, the large space of polynomials can be constructed from a small number 
of representatives via substitutions (see Remarks 12 . 21 13 . 21 and 13 . 5| ) . Theorem 13.101 solves the 
set-theoretic version of Allman's question (up to high numerical accuracy), uses equations 
of lower degree than Friedland's solution [FrilOj . and gives evidence for a conjecture to the 
ideal-theoretical question asked by Allman. 
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One practical interest of the secant variety 04 (P 3 x P 3 x P 3 ) is in phylogenetics, where the 
secant variety is associated to the statistical model for evolution called the mixture model 
of independence models [AR08J. The main motivation to study this particular model is that 
Allman and Rhodes showed in [AR08t Theorem 11] that finding the polynomial invariants 
for this small evolutionary tree would provide all polynomial invariants for the statistical 
model for any binary evolutionary tree with any number of states. 

Note that in this paper we work exclusively over the complex numbers, however, in phylo- 
genetics one is often interested in studying models restricted to the real numbers, the positive 
real numbers, or the probability simplex. Since equations for a given model considered over 
the complex numbers also provide equations for the restricted model, it is natural to start 
with the complex setting and then study the additional necessary equations and inequalities 
imposed by the given restriction. We leave this further study to other works. 

While Allman asks for the generators of the defining ideal of the secant variety, a collection 
of set-theoretic defining equations provides a necessary and sufficient test for membership 
on the model. Very recently Friedland [FrilOj has proved (without a computer) that a 
set of polynomials in degrees 5, 9 and 16 define 04 (P 3 x P 3 x P 3 ) set-theoretically. Indeed, 
Friedland's set of polynomials do (in theory) allow one to test whether a given set of data 
fits the model. Because it uses polynomials in smaller degree, Theorem 13.101 provides a more 
efficient practical membership test for the model. 

On the other hand, Casanellas and Fernandez-Sanchez |CF09] have studied more practical 
issues regarding phylogenetic tree construction using algebraic methods. In particular, they 
point out that for phylogenetic tree reconstruction, the equations coming from the edges of 
the tree (minors of flattenings below) seem to be more relevant than the equations coming 
from vertices (the equations of degrees 5 and 9 are examples of such). 

Our equations in degree 6 are not in the ideal of the equations in degree 5, thus they 
are non-trivial generators in the ideal, and Friedland's result cannot be a set of minimal 
generators of the ideal. We have not found any such obstructions to our result holding 
ideal-theoretically and this leads to a salmon conjecture that the ideal-theoretic version of 
Theorem 13.101 also holds. 

This work was initiated in October 2008 when Bernd Sturmfels asked for a Macaulay2 
readable file of the degree 6 polynomials in the ideal of 04 (P 2 x P 2 x P 3 ). Proposition 
12.11 is a representation theoretic description of these polynomials and corrects minor errors 
in [LM04[ Proposition 6.3], and |LM08[ Remark 5.7]. In Section [2] we give a brief overview 
of how these polynomials were constructed from their representation theoretic description. 
These equations and other ancillary materials for this paper are available in the ancillary 
materials which accompany the arXiv version of this paper or by contacting either author. 

At the December 2008 MSRI workshop on Algebraic Statistics, Oeding presented Conjec- 
ture 13.81 which, when combined with an argument of Landsberg and Manivel, implies our 
main result. This argument is discussed in Section [31 The missing ingredient for the conjec- 
ture was to understand the zero-set of the degree 6 polynomials. Shortly after this workshop, 
the Oeding asked for help from Bates and the Bertini Team. 

The two authors worked together to get the correct mixture of initial input and computing 
strategies in order to find a computation that would finish in a reasonable amount of time. 
Finally on July 12, 2010, a computation that had taken approximately 2 weeks on 8 pro- 
cessors (two 2.66 GHz quad-core Xeon 5410s set up as one head processor and seven worker 
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processors) finished, providing a numerical proof to Conjecture 13.81 Because our calculations 
use numerical approximations, we say that the proof holds up to high numerical accuracy. 
In Section H] we discuss our computational methods and the reliability of this result. 



2. Symmetry and the equations in degree 6 

In this section we recall well-known facts about the variety and equations we are studying. 
The main purpose is to set up notation. The reader who is unfamiliar with these concepts 
may consult [FH91J, or for a more detailed account related to secant varieties see [LM04[ 
ILM081ILW07] or the upcoming [LanlOj . 

Let A, B, C be vector spaces of dimensions a, b, c respectively. The symmetry group of 
a r (FA xPfix PC) is the change of coordinates in each factor GL(A) x GL(B) x GL(C) 
(or when A = B = C there is an additional symmetric group 63 acting and the symmetry 
group is (GL(A) x GL(B) x GL(C)) x ©3 ). Therefore we can use tools from representation 
theory to aid in our search for defining equations. Since much of this work has already been 
done, we only describe the equations relevant for our application. 

The module S d (A* C*D B* ® C*) of degree d homogeneous polynomials on A <g> B ® C has 
an isotypic decomposition (see |LM044 Proposition 4.1]) 



S d {A* ® B* ® C*) = (S ni A* <g> S n2 B* <g> SV 3 C*) ( 



|""l| = |""2 



(-TTl ,7TO ,7TQ 



\n 3 \=d 



where the 7Tj are partitions of d and the multiplicity m 7ri vr2i7r3 is the dimension of the highest 
weight space which can be computed via characters. The modules 

(S Wl A* ® S n2 B* ® S na C*) m *w' 

are called isotypic components, and the individual modules S ni A* <g> S n2 B* (g) S n3 C* are 
irreducible GL(A) x GL(B) x G'L(C)-modules sometimes called Schur modules. 

The ideal of any GL(A) x GL(B) x GL(C)-invariant variety in F(A ® B <g> C) consists 
of a subset of the modules occurring in the isotypic decomposition. If X is a projective 
variety, let X S (X) denote the ideal of homogeneous degree s polynomials in the ideal of X. 
In general, if X is any variety with ideal generated in degree 2 (of which the Segre variety 
is an example), X^a^X)) = for s < k (see [LM04[ Corollary 3.2]), and in particular, 
T s (o-i(FA x FB x PC)) = for s < 4. Also, one can calculate (by checking every irreducible 
module of degree 5 polynomials) that 

X 5 (ct 4 (P 2 x P 2 x P 3 )) = 0. 

In addition, we have found the following. 

Proposition 2.1. Let A = B = C 3 , C = C 4 ; and let M 6 denote the module S 2 , 2 ,2A* 
S 2 ,2,2B* ® S 3tW C*. Then M 6 = X 6 (a 4 (FA xFBx PC)) as GL(A) x GL(B) x GL(C)- 
modules. 

Proof. The module Mq was found by following the ideal membership test described in |LM04j . 
We repeated the procedure outlined in |LM04j as follows. We first decomposed S 6 (A*®B* (g> 
C*) into its isotypic decomposition. This plethysm calculation can be done in one line on 
the program LIE |vLCL92] as 

plethysm([6] ,[1,0,1,0,1,0,0] ,A2A2A3) 
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or by using the procedure "mults" which we implemented in maple and can be found in 
the file "iso_mults.mw," available with our ancillary materials. Next we computed a ba- 
sis of the highest weight space for each isotypic component. We implemented in Maple a 
standard algorithm to compute a basis of the highest weight space in the image of the rel- 
evant Schur functors associated to each module. This implementation is in the file called 
"poly_make_algo.mw" which also may be found with our ancillary materials. A detailed 
exposition of this concept may be found below. We applied this algorithm to each module in 
the isotypic decomposition and we checked to see if any linear subspace of the highest weight 
space of an isotypic component vanished on the variety by direct evaluation. The only module 
which passed this test was Mq, which occurs with multiplicity one in S 6 (A* ® B* cg> C*). □ 

We note that there was some confusion between the statements and proofs in the pre- 
print and the print version of [LM044 Proposition 6.3] as well as in the statement |LM08l 
Remark 5.7], and we believe that Proposition 12.11 corrects this confusion. 

The module 6*2,2, 2C 3 is one-dimensional and as a vector space, the module £3,i,i,iC 4 is 
isomorphic to >S 2 C 4 , which is 10-dimensional. Our construction produces a basis of the 
module M 6 consisting of 10 polynomials which also correspond to the 10 semi-standard 
fillings (strictly increasing in the columns and non-decreasing in the rows) of the tableau 
of shape (3,1,1,1) with the numbers 1,2,3,4. We list these fillings below. The basis of 
polynomials is contained in the file "deg_6_salmon.txt" which is available with our ancillary 
materials as mentioned above. 

Here is a brief overview of an algorithm to construct the polynomials in S ni A* ® S n2 B* g) 
Sn 3 C*. While this algorithm is based on classical methods, we refer the reader to the 
works |Lanl0pOed09aPOed09b] which use similar language for more details. We point out that 
the complexity of any algorithm to compute polynomials from Schur modules will depend 
on dimension and degree. This piece-by-piece algorithm attempts work with the smallest 
dimensional space possible at each step, thus reducing the complexity and increasing the 
chances that the computation will finish in a reasonable amount of time. 

For concreteness, we fix the degree d = 6. The input to the algorithm is the fillings of 
the tableau of shapes 7Ti,7r2,7T3. The first step is to construct a highest weight vector in 
A® G <S> B® 6 ® C® 6 . For this we work one vector space at a time. Suppose a\, 02, 03 is a basis 
of A*. Then a\ ® a\ ® a>i <S> ® a 3 <g> a 3 is a pre-highest weight vector. The Young symmetrizer 

Y m : A* <g> A* <g> A* <g> A* <g> A* <g> A* -> A* <g> A* <g> A* <g> A* <g> A* <g> A* 

is the map that skew symmetrizes the vector spaces A* in positions corresponding to the 
columns of the filling associated to tt\ and then symmetrizes the vector spaces corresponding 
to the rows of the filling associated to t\\. The image of the pre-highest weight vector is a 
highest weight vector of S ni A in (A*) m . We perform the analogous construction in the B* 
and C* factors and take the tensor product of the resulting highest weight vectors. 

The resulting vector we have constructed is in S ni A* (g) S n2 B* <S> S m C*, however it is em- 
bedded in (A*)® 6 ®^*)® 6 ® (C*)® 6 . The final step is to perform the re-ordering isomorphism 
(A*) 06 <g> (S*)® 6 <g> (C*)® 6 -> (A* 5*® C*)® 6 , and then symmetrize the result to arrive at 
a polynomial in S 6 (A* ® B* <g> C*). 
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We computed the 10 polynomials in 62,2,2^* <g) ^2,2, 2B* ® S^i^iC* using the fixed fillings 



_L_ 

5 



1 


4 


2 


5 


3 


6 



for 7Ti and 7r 2 respectively with each of the following fillings for 7r 3 



2 
3 
4 



1 


3 


3 


1 


3 


4 


1 


4 


4 


2 
3 
1 


1 


2 
3 
4 




2 
3 
4 





1 


1 


2 


2 






3 






4 







1 


1 


3 


2 






3 






4 







1 


1 


4 


2 






3 






4 







1 


2 


2 


2 






3 






4 







J_ 

4 



2 3 



J_ 

j2_ 
Ji_ 
4 



2 4 



Notice that up to re-naming the numbers, the fillings for -k^ can be divided into two classes, 
depending on whether the last two numbers in the first row are equal. The four fillings of 
the first class (with the last two numbers in the first row equal) correspond to polynomials 
with 936 terms, whereas the six fillings of the second class correspond to polynomials with 
576 terms. 

Denote by p iJtk , ( 1 < i, j < 3 and 1 < k < 4) a basis of A* <g> B* <g> C* = C 3 ® C 3 ® C 4 . 



takes the polynomial associated to the filling 
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to the polynomial associated to the 



and if m is different than k and I, the swap takes the polynomial associated 
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to the one associated to 
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to the filling LAJ to the one associated to lLI . This additional symmetry could be 
useful for the Bertini computation, however our computation completed without the need 
to implement this symmetry, so we did not use it. We hope to exploit this for future work. 

These fillings produce homogeneous polynomials that are, moreover, homogeneous in 
multi-degree. In general, the multi-degree of a monomial is a collection of vectors [[/^, if], 
[If, I2 , If], [li , 1*2 , 13 , 14]], and is defined on a single variable Xij^ by ifi is (respectively 1) 
for Xij t k ifi^i' (respectively % = i') (If, and 1%, are defined similarly) and the multi-degree is 
defined for monomials by declaring it to be additive over products of variables. For example, 
the following is a sampling of terms in the highest weight polynomial corresponding to the 
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" " " — ^321^113^211^221^134^332 — a;321^1223 ; 231 a; 313^114 + ^211^312^131^1212; 334^223 • • • 

and one finds that this polynomial has multi-degree [[2, 2, 2], [2, 2, 2], [3, 1, 1, 1]]. 

Remark 2.2. Note that when a = b = 3 and c = 4, S 2 ^ 2 A* <g> S 2 ,2,2 B * ® S 3>1X1 C* is 
10-dimensional. When a = b = c = 4, the dimension of ^2,2,2^* ® S 2)2:2 B* <2> S^i^tC* 
increases to 1000, however the basis of this larger space can still be constructed from the two 
polynomials that have 576 and 936 monomials via the type of swap of variables described 
above for the index k in Pij^, but also allowing similar swaps for each of the indices i and j. 

3. Geometric techniques for secant varieties 

Suppose A' C A, B' C B and C C C. Landsberg and Manivel have shown how to take 
equations on ay (PA' x FB' x PC") to equations on ay (PA x FB x PC) and call this procedure 
inheritance |LM04t Proposition 4.4]. 

Subspace varieties contain tensors that can be written using fewer variables. More specif- 
ically, 

Smv j6V (A <8> B ® C) := {[T] G P(A <g> B <g> C) | 3C a ' C A, 

C b ' C 5, C c ' C C, with [T] G P(C a ' <g> C 6 ' ® C c ')}. 

Landsberg and Weyman have shown that Surv,&' iC /(A <g> B <g> C) is normal with rational 
singularities, and the ideal is generated by minors of flattenings |LW07t Theorem 3.1]. Recall 
that a flattening of a 3-tensor in A® 5® C is the choice to view it as a matrix in A® (f?(g> C), 
5 ® (A <g> C) or (A <g> 5) ® C. 

The subspace varieties are important in light of equations because of the fact that 

Sub r)r>r (A (g) _B ® C) D a r (PA x FB x PC), 

and therefore when non-trivial, the ideal of Sub r , r , r gives equations of a r . There is an easy 
test for a module to be in the ideal of a subspace variety, namely S ni A* £g> S n2 B* (g) S n3 C* 
is in the ideal of Sub a /,b', c /(A <S> 5 £g> C) if and only if at least one of the following holds; 
#(7Ti) > a', #(vr 2 ) > 6' or #(^3) > c', where #(•) is the number of parts of the partition. 

Landsberg and Manivel made an important reduction for the salmon problem, which we 
record here. Friedland pointed out that their proof contained an error, which he corrected 
in [FrilOj . Let a, b, c respectively denote the dimensions of A, B, C. 

Theorem 3.1 (Landsberg-Manivel, Friedland). As sets, fora,b,c > 3, 
a 4 (P a_1 x P b_1 x P^ 1 ) is the zero-set of the union of: 

(1) Strassen's commutation conditions, 

M 5 := S(3,1,1)A* <g> 5(2,1,1,1)5* ® £(2,1,1,1) C* 

©£(2,1,1,1) A* ® S (3>ljl) B* <g> £(2,1,1,1^* 

©S(2,1,1,1)A* ® S (2,1,1,1) B * © 'S'O.I.IJC*, 

(2) Equations inherited from a 4 (P 2 x P 2 x P 3 ) ; and 
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(3) Modules in S 5 (A* <S> B* ® C*) containing a f\ 5 , i.e. equations for Sub4 i 4 i 4. 



Note that when a = b = c = 4, the third set of equations is trivial. The key point is that 
we will have a complete description of the set-theoretic defining equations of 04 (P 3 x P 3 x P 3 ) 
as soon as we have the equations of 04 (P 2 x P 2 x P 3 ). 

Remark 3.2. The equations in degrees 5 as well as equations in degree 9 inherited from 
<r 4 (P 2 x P 2 x P 3 ) were found by Strassen |Str83j and were described in terms of certain com- 
mutation conditions. Later, Landsberg and Manivel |LM08] reinterpreted these conditions 
from the geometric and representation theoretic point of view and provided generalizations 
in this language. In |Stu09] one finds a nice description of these equations requiring only 
basic linear algebra. Analogous to our description of the equations in degree 6, here we give 
the representation theoretic description of the polynomials of degree 5. 

Note also that when a = b = c = 4, M5 is a 172 8- dimensional irreducible G-module, 
for G = (GL(4) x GL(4) x GL(4)) x & } . A natural basis of M 5 can be constructed as in 
the previous section. For this we need to give the fillings for the triple of Young diagrams 
corresponding to the partitions (2, 1, 1, 1), (2, 1, 1, 1), (3, 1, 1). We note that up to permu- 
tation, there is just one equivalence class for the fillings of the diagram for (2, 1, 1, 1) with 
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There are three equivalence classes for the fillings of the diagram for 
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> £(2,1,1,1) -B* 



1 5(3,1,1) C* 



WC 



(3, 1, 1) with representatives 

Therefore to construct representatives for a basis of £(2,1,1,1)^* 
fix the representative filling for (2, 1, 1, 1) in both instances, and we let the filling for (3, 1, 1) 
vary over the three representatives. Thus we construct three polynomials, one for each 
representative filling of the diagram for (3, 1, 1) and respectively, these polynomials have 180, 
360 and 540 monomials. A basis of polynomials for one of the 3 isomorphic modules in M5 
is contained in the file "deg_5_salmon.txt" with our ancillary materials. After constructing 
these three polynomials, the rest of the polynomials in the basis of M5 can be constructed 
by the substitutions and swaps of variables as mentioned above (see the discussion above 
Remark 

Another important result for the salmon problem is from Strassen, which has been rein- 
terpreted in representation theoretic language in |LM 08j. 



Theorem 3.3 ( |Str83j ). The ideal of the hyper surf ace <r 4 (P 2 x P 2 x P 2 ) c P 26 is generated 
in degree 9 by a nonzero vector in the 1- dimensional module 

S(3,3,3)C 3 ® %3,3)<C 3 <g> 5 (3 ,3,3)C 3 . 

Let Mg denote the inherited module £(3.3,3)C 3 Cg) £(3,3,3) C 3 ® £(3,3,3) C 4 . Inheritance implies 
that M 9 e Z(a 4 (P 2 x P 2 x P 3 )). 

Remark 3.4. Suppose [T] £ ¥(A <g> B ® C), with dim(A) = 3. Then write T = a 1 ®T 1 + a 2 ® 
T 2 + 03 Cg> T 3 , where the T, are b x c matrices in B <g> C and the basis of A. 
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Strassen described his equation in degree 9 as follows. On an open set one may assume 
that T\ is invertible. Then consider the polynomial 



det^) 2 det(T 2 Tr 1 T 3 - T 3 T^T 2 ) 



Strassen showed that this polynomial is irreducible, of degree 9, and vanishes on 04 (PA x 
FB x PC). 

A useful reformulation by Ottaviani of Strassen's equation is the following (see |LQ1CH 
IOtt07] ). As before write T = a\ <g> T\ + a 2 ® T 2 + a 3 ® T 3 . Here one does not require any of 
the slices T 1; T 2 , T 3 to be invertible. Construct the block matrix 

T 3 -T 2 
= ( -T 3 Tx 
T 2 -Ti 

One checks that ip T is linear in T, and that if [T] e Seg(FA x FB x PC) then Rank(i/j T ) = 2. 
Therefore if [T] is a general point in <7k{FA x FB x PC) it can be written as the sum of 
k points on Seg(FA x FB x PC) so Ranker) < 2A; by the sub-additivity of matrix rank. 
In particular, in the case dim(A) = dim(S) = dim(C) = 3, the 9x9 determinant det(?/>r) 
gives a non-trivial equation for 04 (PA x FB x PC), which is also Strassen's equation. This 
polynomial has 9,216 monomials. Note that ipT is not a skew-symmetric matrix unless the 
matrices T{ are symmetric, otherwise any odd-sized determinant would vanish identically. 

Remark 3.5. In the case that a = b = 3 and c = 4, as a vector space, Mg is isomorphic to 
5 3 C 4 so dim(M 9 ) = 20. When the highest weight vector of a module has a determinantal 
representation (as in the case of Mg), it is typically much faster to compute a basis of the 
module from the highest weight vector using lowering operators. (Lowering operators are 
standard in the theory of Lie algebras, but are not the focus of this work. We refer the 
interested reader to |Oed08t Section 3.4] for an explicit treatment of this method.) Using 
this method, we found that the natural basis of Mg consists of polynomials with 9, 216 or 
25,488 or 43,668 monomials. This basis is a 23Mb text file of polynomials, too large to 
include with our ancillary files due to the restrictions of the arXiv, but may be obtained 
from either author. As in Remarks 12.21 and 13. 2[ these polynomials can be associated to 
representative polynomials, depending on fillings. In the A and S-factors, the diagram for 



(3,3,3) can only have one semi-standard filling, namely 



are three classes of fillings, namely 
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In the C-factor, there 



These fillings yield the 

representative polynomials consisting of 9,216 or 25,488 or 43,668 monomials respectively. 
The rest of the polynomials in a basis of Mg can be constructed by the substitutions and 
swaps described in our treatment of M 6 (see the discussion above Remark |2~2]) . 

Alternately, a basis of Mg can be constructed via Ottaviani's formulation. They are derived 
from the condition that the now 9 x 12 matrix appearing in (CQ) have rank 8 or less. However, 
the space of 9 x 9 minors of ipT is no longer irreducible when a = b = 3 and c = 4. Namely 
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the space of 9 x 9 minors of the 9 x 12 matrix ipT is the following representation 

63,3,3^* © 63,3,3^* <g> 6 3 ,3, 3 C* 
©64,3,2^* © 33,3,35* ® 6 3 , 3 ,2,lC* 
©6 5 ,2, 2 A* © 6 3 ,3, 3 P* © 6 3 ,2,2, 2 C*. 

There are three equivalence classes of maximal minors of ipT depending only on the column 
index / of the maximal minor of Aj(i/jt)- Let P = (Pi, P2, P3) be the partition of the set 
{1, . . . , 12} into three sets Pi = {1,2,3,4}, P 2 = {5,6,7,8}, P 3 = {9,10,11,12}. The 
representation 6*3,3,3 A* © S^^B* cg> 63,3,3c* is associated to the minors Aj(t/V) such that 
\I fl Pi\ =3 for i = 1, 2, 3. This condition precisely forces the minor of ipT to be constructed 
with 3x3 submatrices of Ti,T 2 and T3. The representation 5 4 ,3 )2 A* ® S 3 ^ 3t3 B* <S> 63,3,2,1c* 
is associated to the minors Aj(^r) such that \I fl Pi| =4, \I D P 2 | = 3, |J D P3I = 2. The 
representation 65, 2,2 A* © 83,3,38* © 5*3, 2,2,2c* is associated to the minors A/(^ T ) such that 

|/nPi| =4, |/nP 2 | =4, |/nP 3 | = 1.' ' 

Note that the symmetry implied by the fact that A and B have the same dimension allows 
us to reverse the roles of A and B to find two more modules in the ideal, namely the two 
modules S 3 ^ 3 A* <g> 6 4 , 3 , 2 P* <g> 6 3 , 3 ,2,iC* and £3,3,3^* <g> 6 5 , 2 ,2P* © 63,2,2,2c* must also vanish 
on a A {FA x PP x PC).' 

While we have described five modules of degree 9 equations which vanish on 04 (PA x 
PP x PC), we only use the module M 9 = £3,3,3 A* <g> 6 3 , 3 , 3 P* <g> 63,3,3c* along with M 6 
described above for our set-theoretic defining equations. We can conclude that (M 9 ) <£_ (M 6 ) 
by analyzing the shapes of the partitions involved. More specifically, in the C-factor the 
partition (3,3,3) only has 3 parts, but if S ni A* ® S^B* ® S n3 C* is a module in the ideal 
generated by M 6 then ir 3 must have at least 4 parts. However this argument fails for the 
other four degree 9 modules so it is possible that these equations are in the ideal generated 
by Mq. Moreover our set-theoretic result implies that it must be the case that the other 
degree 9 modules are in the ideal generated by M 6 (up to high numerical accuracy). 

Example 3.6 ( |Fril0j ). Friedland has shown that the known equations in degree 9 are not 
sufficient to define <7 4 (PA x PP x PC) set-theoretically when dim(A) > 3, dim(P) > 3 and 
dim(C) > 4. We thank J.M. Landsberg for the following clarification of Friedland's example. 
Consider the point 

P = (ai®6i + a2©62)©Ci + (ai©fci + a2©^3)©C2 + (ai©fti + a3©^2)©C3 + (ai©6i + a 3 (8)63)(8)C4. 

The span of {01,02,03} C A and the span of {61,62,63} C P are both no more than 3- 
dimensional, so P is a zero of M 5 since the representations S ni A* <g) S^B* <g) S n3 C* in M 5 
each have either | vri | = 4 or =4, and therefore the respective Schur functor S 7Ti with 
|7Tj| =4 will annihilate a 3-dimensional subspace. One finds that ^t(P) has rank 8 and 
therefore P is a zero of Mg. However P is not a point of a^(FA x PP x PC). This geometric 
argument implies that more polynomials are needed than just the degree 5 and 9 equations. 
For this, Friedland produces equations of degree 16 which do not vanish on P. On the other 
hand, P is not in the zero set of M 6 , so M 6 is sufficient to rule out the possibility of points 
of the same form as P to have border rank 4. Therefore, one could repeat Friedland's proof, 
modifying the argument where he uses degree 16 equations with these degree 6 equations 
and thus obtain a new result, and a computer-free proof of Theorem 13. 101 
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Remark 3.7. To construct a basis of the 8000 dimensional space S^^C 4 ($ S^^^C 4 (g) 
5(3 3 3 ) C 4 , one can repeat the lowering operator procedure. Since these polynomials are very 
complicated, our experience is that, in practice, one should use the degree 9 equations in their 
determinantal form. In particular, to check if a point z vanishes on all of the polynomials in 
>S'(3,3,3) ( C 4 g> S^^^C 4 g) S^^C 4 , it is more efficient to first construct the matrix in ([T|) for 
the point z and check that the determinant vanishes. Then repeat this test for all allowable 
changes of coordinates, in other words, for every g G GL(4) x GL(4) x GL(4) construct 
the matrix in ([[]) for g.z and check that the determinant still vanishes. (This is sufficient 
because our module is the span of the orbit of a single polynomial.) Moreover, if one only 
wants a quick check that z is in the zero-set with high probability, it suffices to check that 
g.z is in the zero-set for a random g. In this quick test, a non- vanishing result is certain, but 
vanishing must be re- verified with an exact (non-randomized) test. 

Since (2,2,2) has 3 parts, and (3,1,1,1) has 4 parts, M§ must vanish on the subspace 
varieties Sub^^ U Sub 3j2j 4 U Subs^^. Also, note that two of these subspace varieties are al- 
ready contained in the secant variety, namely 04 (P 2 x P 2 x P 3 ) d Sub^^ U Suba^. Indeed, 
if x G Sub 2 ,3,4, there exists A' C A such that dim (A') = 2 and x G F(A' <g> B <g> C). But in 
this case F(A' <g> B ® C) = a A (FA' x PS x PC) C a A (FA x PS x PC). The same argument 
is repeated for Sub 3i2 ,4. 

If M is a set of polynomials, let V(M) denote the zero-set of M. Based on the above 
evidence, we make the conjecture 

Conjecture 3.8. As sets, 

V(S (2>2 , 2) C 3 ® 5( 2 , 2 , 2) C 3 <g> S (3Xltl) C A ) = a A (P 2 x P 2 x P 3 ) U Sub 3 ,3 i3 . 

Computation I4TT1 below verifies that Conjecture 13.81 is true up to high numerical accuracy. 

Theorem 3.9 (Corollary to Computation KH) ■ Let A = C 3 ,B = C 3 , C = C 4 . Up to high 
numerical accuracy, the secant variety 04 (FA x FB x PC) is defined set-theoretically by 

M 6 = S (2 ,2,2)A* <g> S (2 ,2,2)B* ® S {3ilAA) C* 

Mg = 5(3,3,3)^* ® S^ 3 ,3)B* (g) S(^ 3 j)C*. 

Proof. By Proposition 12.11 and by Strassen's Theorem 13.31 combined with inheritance we 
know that both M 6 and M 9 are in the ideal of a 4 (FA* x FB* x PC*). So we know that 
a 4 (FA x PS x PC) C V(M 6 © M 9 ). 

For the other inclusion, select a point z in the common zero locus of Mq and Mg. Since 
z G V(Me), Conjecture 13.81 says that either z is on the secant variety, in which case we are 
done, or z is on the subspace variety. In the latter case, let C'cCbea 3-dimensional vector 
space so that z G F(A®B®C). Then z is a zero of M 9 = S {3)3)3) A*(g)S {3A3) B*(g)S {3A3) C*, and 
therefore is also a zero of the polynomials in the restriction S( 3>3>3 )A* <g> S^^B* (g> ^^sjC*. 
So by Strassen's Theorem 13. 3[ 

2 G (T 4 (PA x PS x PC) ^ a 4 (P 2 x P 2 x P 2 ). 

We conclude because we have the obvious inclusion 

a 4 (FA x PS x PC) C a A (FA x PS x PC). 

□ 
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We used numerical algebraic geometry, specifically Bertini, to compute the decomposition 
of the zero set V(Mq) into irreducible varieties. We outline this computation in the next 
section. However, if one were to prove Conjecture 13. 8[ then the qualifier "with high numerical 
accuracy" may be removed from the statement of Theorem 13.91 

Recall that the Landsberg-Manivel-Friedland Theorem 13.11 above said that set-theoretic 
defining equations of a 4 (P"" 1 x P^" 1 x P c ~ 1 ) with a,b,c > 3 will be known as soon as 
set-theoretic defining equations of cr 4 (P 2 x P 2 x P 3 ) are known, and this is the content of 
Theorem 13.91 Therefore we can restate the immediate consequence of combining Theorem 
13.11 with our computations: 

Theorem 3.10. As sets, for a, b, c > 3, up to high numerical accuracy, 0-4 (P a_1 x P 6 ^ 1 x P^ 1 ), 
is the zero-set of: 

(1) Strassen's commutation conditions, 

M 5 := Sp A)1) A* <8> 5 (2 ,i,i,i)5* ® 5(2,1,1,1) C* 

©S(2,1,1,1)A* <g> 5(3,1,1)5* ® 5 (2 ,1,1,1)C* 

®S {2 ,i,i,i)A* (8) 5(2,1,1,1)5* <g> 5 (3 ,i,i)C*, 

(2) equations inherited from cr 4 (P 2 x P 2 x P 3 ), 

Mq = 5( 2 ,2,2)^4* <S> 5(2,2,2)5* <g) 5( 3 ,i,i,i)C* 
M 9 = 5(3,3,3)^4* ® 5(3,3,3)5* ® 5(3,3,3) C*, 

(3) and modules in S 5 (A* £g> 5* <g> C*) containing a /\ 5 , i.e. equations for Sub4,4,4. 

Remark 3.11. The qualifier "up to high numerical accuracy" can be removed if one uses 
Friedland's argument [FrilOj modified by our computations, as mentioned in Example 13.61 

4. Results using numerical algebraic geometry 

In this section, we provide a brief overview of the basic methods of numerical algebraic 
geometry; references for further details are provided. We then describe the results of the run 
establishing the main result of this article and conclude with a short discussion regarding the 
reliability of numerical algebraic geometry methods and, more to the point, the reliability of 
this result. 

4.1. Brief overview of numerical algebraic geometry methods. Given generators of 
an ideal of C[x\, . . . ,xn], the methods of numerical algebraic geometry will produce a nu- 
merical irreducible decomposition for the associated variety X C C^. In particular, for 
each irreducible component Z of X, these methods will produce deg Z numerical approxi- 
mations (to any number of digits) of generic points on Z. The end result is a catalog of all 
irreducible components of X, each indicated by a set of witness points on the component 
(together referred to as a witness set for the component), its dimension, and its degree. 

The core method of numerical algebraic geometry is homotopy continuation, a method for 
approximating the complex zero- dimensional solution set of a polynomial system. The basic 
idea of homotopy continuation is to cast the given polynomial system F as a member of 
a parameterized family of polynomial systems, one of which, G, has known solutions or is 
otherwise easily solved. If done correctly, the solutions of G will vary continuously to those 
of F as the parameters are varied appropriately. By tracking these paths numerically (using 
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predictor-corrector methods), one will arrive at numerical approximations of all complex 
zero-dimensional solutions of F. There have been many technical advances in this area 
that contribute heavily to the reliability of these methods. See |SW05tlLi03j for general 
references and [BHSW 08, BH SW09] regarding the use of adaptive precision methods for 
added reliability. 

Pairing homotopy continuation with the use of hyperplane sections, monodromy, and a 
few other methods described fully in |SW05] yields the numerical irreducible decomposition. 
Briefly, a d- dimensional irreducible algebraic set in will intersect a generic codimension 
d linear space in a set of points. This statement about genericity (along with similar as- 
sumptions of genericity throughout numerical algebraic geometry) is the reason for referring 
to these methods as probability-one methods, as described further below. 

The computation of a numerical irreducible decomposition begins by searching for codi- 
mension one irreducible components (by adding N — 1 linear polynomials to the set of gener- 
ators and solving for zero-dimensional components via homotopy continuation), followed by 
codimension two components, etc. Once this sweep through all possible dimensions has been 
completed, we have a superset of the desired numerical irreducible decomposition, since a lin- 
ear variety of codimension d will intersect any component of dimension d or higher. Sommese, 
Verschelde, and Wampler (and others) have developed methods for removing points in the 
"wrong dimension," i.e., those discovered while searching for components in dimension d 
which actually lie on higher- dimensional components, called junk points. They have also 
developed algorithms for performing pure-dimensional decompositions to yield witness sets 
on each irreducible component (instead of the initially- found witness sets for the union of 
all equidimensional irreducible components). See |SW05j for further details. 

There are three main software packages in this field: Bertini |BHSW10bj . HOM4PS- 
2.0 |LLT10j . and PHCpack |Verl0j . Each package has various benefits over the others |BHSW10a] . 
Since Bertini is typically the most efficient package for large, parallel, positive-dimensional 
problems as well as the package with the most reliability and precision features, we used 
Bertini in our computations for this article. In fairness, it should also be noted that Bates 
is a Bertini developer. 

4.2. Numerical results for the Salmon Problem. 

Computation 4.1. Up to 10 digits of accuracy, the zero-set of the 10 polynomials in a basis 
of Mq (defined above) has precisely two irreducible components. One, in dimension 31, has 
degree 345. The other, in dimension 29, has degree 84- 

Indeed, cr 4 (P 2 x P 2 x P 3 ) is non-defective and has dimension 31 |AOP09t Theorem 4.6]. 
It is also straightforward to check that Subs^^ has dimension 29, and by the pigeon-hole 
principle, these must be our components in the zero-set of M 6 . Though these dimensions are 
sufficient information to identify our varieties, as additional information, we find that this 
secant variety has degree 345 and the subspace variety has degree 84. 

Proof. The conclusion comes from the results of a calculation on Bertini |BHSW10b] using 
approximately 2 weeks of computing time on 8 processors, using tight controls including small 
tracking and final tolerances (10~ 10 or smaller), adaptive precision numerical methods, and 
a variety of checks and error controls built into Bertini (such as checking at t = 0.1 that no 
paths have crossed). The output of our computation is included in the files "main_data.txt" 
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and "screen_out.txt" , which may be obtained with the other ancillary materials as mentioned 
above. □ 



Remark 4.2. After the submission of the first draft of this paper, Friedland provided a 
counterexample to |LM08t Proposition 5.4], thus invalidating the proof of [LM08t Corol- 
lary 5.6]. Because we quote this result in the present paper, we tried to use numerical 
methods to find out what could be true. In particular, we tried to find a corrected statement 
for |LM08t Proposition 5.4] by computing the zero-set of S^i^A* <g) S(2,i,i,i)5* <S> £(2,1,1,1) C* 
using Bertini. This proved to be a very expensive computation. This module, in its smallest 
form, has a basis of 96 degree 5 polynomials in 48 variables (see the file "deg_5_salmon.txt" 
in our ancillary files). After one month of computational time on 72 processors, we had only 
completed the first 10 codimensions, tracking up to 2 million paths for each codimension. An 
easy geometric argument implies that the smallest possible component has projective dimen- 
sion 8, indicating that we were very far from completing the computation. In the mean time, 
a second version of [FrilOj appeared, with a corrected proof for [LM08[ Corollary 5.6]. Due 
to limited computational resources and time, we decided to abandon further computation 
and accept the computer- free proof of [LM081 Corollary 5.6] in |Fril0j . 

4.3. Reliability of this result. Can Computation I4TT1 be accepted as absolute proof? No, 
unfortunately, it may not. However, this numerical computation gives extremely strong 
evidence that Computation 14.11 is indeed true even without the phrase "with high numerical 

accuracy" . 

There are two types of approximations that are used in order to compute the numerical 
decomposition of a zero-set. One is the choice of a set of random hyperplanes which cut 
the space and allow one to look for zero-dimensional solutions to a set of equations. The 
other type of approximation is the numerical homotopy continuation method which actually 
searches for the zero-dimensional solutions. 

The choice of random hyperplanes amounts to the choice of random numbers from a Zariski 
open, dense set S of some parameter space rather than choosing some set of points in the 
complement of S. Since the complement of S is an algebraic set, we know that it must have 
positive codimension, making it a set of measure zero for any reasonable choice of measure. 
Thus, the set of hyperplanes that fail in that they would cause us to miss a component in the 
zero-set has measure zero, and we say that the choice of hyperplanes will yield the correct 
result with probability one. 

The second type of approximation that is done in this type of computation is the heart of 
Bertini and is thoroughly described in [SW05] . Bertini allows one to set desired accuracy to 
arbitrary levels, and any computational errors (such as path crossing) are reported. Further, 
Bertini has additional features such as adaptive precision path tracking, which increase 
security, |BHSW08llBHSW09] . 

The run for this article used a special equation-by-equation algorithm called regenera- 
tion |HSWarj . The run required the following of more than 200,000 paths and there were 
no path failures and no crossed paths detected. In addition, there were no errors in the 
monodromy or trace test procedures. The numerical output of our run is contained in the 
files "main_data.txt" and "screen_out.txt" with our ancillary materials. 

We cannot conclude with unquestionable certainty that Computation 14.11 holds uncon- 
ditionally, but we can state with an extremely high level of confidence that it is correct. 
Motivated by this result, we hope to find a direct argument to prove Conjecture 13.81 
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4.4. Numerical vs. symbolic computation. Finally, one might wonder why we chose 
to use numerical methods to test this conjecture rather than symbolic methods that will 
provide certainty. The main reasons are simple: time and space. Regarding time we expect 
that without additional ideas to reduce the difficulty of computation, a related calculation 
using symbolic methods should take at least eight times as long as the calculation in Bertini 
because Grobner basis algorithms are not completely parallelizable (but for an example of 
recent progress on this front see [Kre09j ) . In fact, based on the timings from an ongoing 
benchmarking project between the Bertini and Singular jDGPSlO] development teams, we 
suspect that any symbolic computation will actually take far more than eight times as long. 
Regarding the issue of space we must consider data storage at intermediate stages. While the 
initial input and final result may be relatively small, Grobner basis algorithms typically must 
store large intermediate results for subsequent calculations. On the other hand, homotopy 
continuation algorithms require a trivial amount of extra data in intermediate stages. Indeed, 
the amount of memory used grows linearly with the number of paths tracked (simply because 
the final point on each path must be stored). Bertini is thus much less likely to fail due to 
memory constraints. 

Finally, one could also hope for a (symbolic) certificate of the validity of results obtained 
by numerical methods. At the EACA School in Tenerife, Spain, Wolfram Decker told us 
that the development of such certificates is among the current goals of the Singular team, 
and we hope to be able to use this feature in future work. 
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